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We derive the molecular properties for a sample of 301 Galactic H ll regions 
including 123 ultra compact (UC), 105 compact, and 73 diffuse nebulae. We 
analyze all sources within the BU-FCRAO Galactic Ring Survey (GRS) of ^"^CO 
emission known to be Hll regions based upon the presence of radio continuum 
^ . and cm-wavelength radio recombination line emission. Unlike all previous large 
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area coverage ^^CO surveys, the GRS is fully sampled in angle and yet covers 
~ 75 square degrees of the Inner Galaxy. The angular resolution of the GRS 
^ , (46") allows us to associate molecular gas with Hll regions without ambiguity 

and to investigate the physical properties of this molecular gas. We find clear 
CO/Hll morphological associations in position and velocity for ~ 80% of the 
, nebular sample. Compact H ll region molecular gas clouds are on average larger 

OO '. than UC clouds: 2'.2 compared to 1'.7. Compact and UC Hll regions have very 

similar molecular properties, with ~ 5 K line intensities and ~ 4kms~^ line 
widths. The diffuse Hll region molecular gas has lower line intensities, ~ 3 K, 
and smaller line widths, ~ 3.5kms~^. These latter characteristics are similar 



OO . 

O . to those found for quiescent molecular clouds in the GRS. Our sample nebulae 

thus show evidence for an evolutionary sequence wherein small, dense molecular 
^ , gas clumps associated with UC H ll regions grow into older compact nebulae and 

^ ! finally fragment and dissipate into large, diffuse nebulae. 

Subject headings: Hll regions — ISM: abundances, clouds, atoms, evolution, 
lines, and bands, structure — radio lines: ISM 
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INTRODUCTION 



In the classical view of Hll region formation OB stars form inside Giant Molecular 
Clouds (GMCs). Any newly formed OB star emits extreme ultraviolet (EUV) radiation 
(> 13.6 eV) and ionizes the surrounding medium of a molecular cloud, creating an H ll region. 
The ionizing photons have more energy than is necessary to ionize the gas, however, and 
thus some excess energy heats the ambient gas. Because of the pressure difference between 
the cold, natal molecular cloud (~30 K) and the Hll region (~10'^K), the ionization front 
expands into the molecular cloud. 

A photodissociation region (PDR) exists beyond the ionization front. The PDR is 
a zone where photons of energies lower than 13.6 eV ionize elemen ts with low ionization 
potentials and dissociate molecules (see iHoUenbach &: Tielens 1119971 . for a review of dense 
PDRs). Within the PDR there are three important boundaries. At Ay — 1, a dissociation 
front exists where H2 replaces H as the dominant species. Further from the emitting star, at 
Ay — 4, there is a boundary between C"'"/C/CO where carbon becomes stored predominantly 
in the CO molecule. The final boundary of the PDR occurs at the O/O2 transition. Ay — 10. 
The PDR is therefore a transition region between the H ll region and the molecular cloud; it 
is completely ionized at one boundary and completely molecular at the other. Since all Hll 
regions should produce PDRs, in this simplified view all H ll regions should have associated 
molecular gas when they first form. As the nebulae evolve, however, the OB stars can travel 
far enough to leave the environs of their natal molecular clouds. 

Complicating this simple scenario is the fact that molecular clouds are clunipy and inho- 



moge neous on all scales (e.g.. iFalgarone &: Phillips Ill996l : iKramer et al. Ill998l : ISimon et al 



20011 ). In a clumpy medium, EUV photons can penetrate to different depths, creating non- 
spherical ionization fronts. Additionally, H ll regions evolve by moving away from their natal 
cloud environments and displacing the local gas. Understanding the complicated interaction 
between young stars, H ll regions, and molecular gas is crucial to the study of massive star 
formation and the impact massive stars have on their environment. 

The molecular component of H ll regions has been studied in detail by many authors. 



Most studies, however, ha ve focused mainly on ultra-compact (UC) H ll regions (e.g.. lChurchwell. Walmsley. 
I990I : iKim &: Koo I l2003l ) . When compact H ll regions were observed in CO, they were of- 



ten only observed with single pointing s using single-dish telescop es (e.g.. lBrand et al. Ill984 



Whiteoak. Otrupcek. &: Rennie Ill982l : iRusseil fc Castets II2004I ). As our results show, CO 
gas is frequently offset from the nominal position of an Hll region. A completely sampled 
map is required to understand fully the dynamics and properties of molecular gas associ- 
ated with H II regions. Furthermore, a sample with H ll regions at all evolutionary stages is 
necessary to understand how this interaction progresses as the H ll region ages. 
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Here we describe a large scale study of the molecular properties of diffuse, compact, and 
ultra-compact Galactic H ll regions using fully sampled CO maps. We trace these molecular 
properties using the ^^CO J = 1 —>■ emission mapped by the Boston University-Five 
College Radio Astron omy Observatory (hereafter BU-FCRAO) Galactic Ring Survey (GRS: 



Jackson et al. II2006I ). CO is an excellent tracer of molecular material because of its high 



abundance and the fact that it is rotationally excited at densities common in molecular 
clouds f>500cm~^). 



2. Milky Way Surveys 

2.1. Galactic H II Regions 

Our Galactic H ll region source sample is compiled from the H ll region radio recombi- 
nation line (RRL) catalog of iLockman I (119891 ) (hereafter L89), which consists of nearly 500 
RRL observations with a 3' beam at positions of known continuum emis sion. These con 
tinuum emission sources were drawn from t he 5 GHz continuum survey of lAltenhoff et al. 
( 1l979l ). All positions in lAltenhoff et al. I (119791 ) with a peak flux density greater than 1 
Jybeam"^ and many sources with peak flux densities down to 0.8 Jybeam"^ were observed, 
unless the source was i n a ver y confused region or was a known supernova remnant (SNR). 
Since lAltenhoff et al. I Jl979h has coverage from / = 357?5 to 60° and |6| < 1°, the L89 
catalog is complete down to a flux limit of at least 1 Jy at 5 GHz in this region of the 
Galaxy. This sky coverage entirely overlaps and is much larger than the GRS survey zone 
(discussed in §2.2p . so the L89 catalog is complete at least down to IJybeam"^ over the 
extent of the GRS. The sources in L89 are the classical Hll regions first detected in the 
1970's at cm-wavelengths using RRLs. Although they are often misidentified, these are not 
the ultra compact, high density Hll regions that are being studied mostly using radio inter- 
ferometers. Below, to distinguish these classical Hll regions from UC nebulae, we will call 
the L89 nebulae "compact" Hll regions. 



Our nebular sample also includes diffuse H ll regions from lLockman. Pisano. fc Howard 

( 119961 ) (hereafter L96). This catalog consists of RRL measurements at 6 cm (~ 6' beam) and 
9 cm (~ 9' beam) toward 13 faint , extended continuum sources. Most of these sources are 



drawn from lAltenhoff et al. I (119791 ) and have peak flux densities greater than 0.5 Jy beam ^. 
The L89 survey is actually the pilot study of this diffuse sample; it contains 40 diffuse sources. 

Recently a large catalog of 1442 Galactic Hll re gions was compiled from 24 published 
studies of Galactic Hll regions (jPaladini et al. II2003I ). These Hll regions were found using 
single-dish, medium resolution (few arcminute beamwidths) observations. For our purposes. 
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however, this catalog is not useful because nebular positions are given to an accuracy of only 
6'. 



Fin ally, our sour c e sample contains UC Hll regio ns taken from IWood &: Churchwell 



f 1989a l. lKurtz et al. I (ll994[ ). IWatson et al. I(l2003f ). and lSewilo. et al. Ifl2004l ). The lWood fc Churchwell 
(ll989al ) UC sources were selected by: (1) the presence of a small or unresolved radio 
source, (2) a spectru m consistent with f ree-free emissio n , and (3) strong FIR ernission . 



The UC sources from Kurtz et al. I (Il994[). IWatson et al I fl2003h . and lSewilo. et al. I (120041 ) 



were selected based on the IWood fc Churchwell I (1l989bl ) criteria: IRAS flu x density ratios 



l oR(Fo jF^y) > 0.57, log(F6o/Fi2) > 1.30, and Fioo > 1000 Jy (> 700 Jy for IWatson et al. 
( 120031 )). where Fx is the IRAS flux density at A /xm. 



2.2. The "CO Galactic Ring Survey 



We use the BU-FCRAO ^^CO Galactic Ring Survey datc0 Jjackson et al. I bood ) to 
characterize the molecular properti es of all H ll regions in th e GRS. The GRS traces th e 
5 kpc molecular ring discovered by iBurton et al. I (119751 ) and IScoville fc Solomon I (119751 ). 
This annulus of enhanced CO emission dominates the inner Galaxy's structure and harbors 
most of the Galaxy's star formation regions. The GRS sky coverage spans 18° < / < 55° and 
\b\ < 1°. Additional, incomplete sky coverage is available for 14° < / < 18° over the same 
latitude range. The GRS covers a total of 74 square degrees. The GRS maps the distribution 
of emission from the J = 1 — (z/q = 110.2 GHz) rotational transition of ^^CO. The ^^CO 
isotopologue is ~ 50 times less abundant than ^^CO and hence has a much smaller optical 
depth. This decreased optical depth yields smaller line widths and gives a cleaner separation 
of individual velocity components along any specific line of sight compared to previous ^^CO 
surveys. 

The GRS, with a spectral resolution of 0.21 km s^^, an angular resolution of 46", and a 
22 " angular sampling, improves upon all previous large scale CO surveys. It is the only fully 
sampled (in solid a ngle) large scale ^^CO survey extant. The GRS improves upon the Bell 



Labs ^^CO survey (ILee et al. II2OOII ) that has a spectral resolution of 0.68 km s ^, an angular 



resolution of 103" and 180" angular sampling. The GRS also has better resolution than 
previous ^^CO survey s. For example, the University of Massachusetts Stony Brook survey 
( Sanders et al. Ill986l ) has a spectral resolution of 1.0 km s~^, an angular resolution of 45" and 
180" angular sampling. The Columbia/CfA ^^CO survey (jPame. Hartmann &: Thaddeus 
2OOII ) has a spectral resolution of 0.18kms^^, an angular resolution of 450" and 225" to 



Data available at 



littp://www. bu.edu/galacticring/ 
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450" angular sampling. Compared to GRS, all of these surveys are severely undersampled in 
angle. The GRS maps have spectra observed at positions separated by ~ | the telescope's 
beamwidth (0.48 HPBW, actually). This, together with high spectral resolution allows us 
to separate individual ^^CO components cleanly. Thus we can for the first time study the 
molecular properties of Galactic H ll regions free from angular sampling bias. 



3. H II Region Source Sample 



Here we study the molecular properties of all known Hll regions in the zone mapped 
by the GRS. Our sample of 301 Hll regions contains 123 UC, 105 compact, and 73 diffuse 
nebulae. For UC nebulae without RRL measurements in the original papers, we compile 



RRL velocities from lAfflerbach et al. I (119961 ) and lAraya et al. I (120021 ) or, if the UC and 
compact sources were co-spatial, from L89. The compact nebulae are from the L89 catalog. 
The majority of our diffuse regions are from L96, but a small number are from the pilot 
survey of diffuse regions in L89. This final nebular sample results from further vetting of 
the §2 sources using data from several recently completed Galactic scale sky surveys that 
overlap the GRS zone. 

We verify the existence, classification, and position of each Hll region by examining 
the radio continuum and infrared emission at its nominal position. For th e radio continuum 
emission, we primarily use the 21cm VLA Galactic Plane survey (VGPS: IStil et al. II2006I ). 
In addition to the 21 cm Hi line emission data cubes, the VGPS generated 21cm continuum 
images over the range 18° < / < 67°, |6| < 1 — 2 with ~ 1' resolution. In addition to the 
VLA measurements, the VGPS used a Green Bank Telescope 21 cm survey to provide the 
zero spacing data. The VGPS is therefore sensitive to both large and small scale emission. 
We also use the 20 cm data from the Multi- Array Galactic Plane Imaging Survey (MAGPIS: 



Helfand et al. II2006I ). These data were collected with the VLA operating in B-, C-, and 
D-configurations and have a resolution of ~ 6". The 20 cm MAGPIS data cover the range 
5° < / < 48?5, |6| < 0?8. We find MAGPIS to be best suited for verifying UC nebulae, 
while the VGPS is better suited for verifying compact and diffuse Hll regions. For the 
infrared emission, we use the 8 nm data from the Ga lactic Legacy Infrared Mid-Plane Survey 
Extraordinaire (GLIMPSE: iBenjamin et al. 1 12003| ) and 24 /zm data from the MIPS Inner 
Galactic Plane Survey (MIPSGAL: Carey et al. 2008 in preperation) . Both infrared surveys 
have coverage beyond the extent of the GRS. 

For inclusion in our sample, we require a continuum peak at the position of each Hll 
region. The UC regions were identified by their IRAS colors and the compact and diffuse 
regions were located in radio continuum maps with 2'.Q resolution. Both of these identification 
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methods have some level of error that can be reduced through correlation with high resolution 
radio continuum data. For UC regions, continuum observations are necessary to confirm 
that the nebula is an Hll re gion, and not a de n se pr otostellar clump. The sources from 



Wood fc Churchwell I (Il989al ) and iKurtz et al. I (Il994j ) were confir med to be UC r egions 



wit h VLA continuum ob servations, but the majority of sources from IWatson et al. I (120031 ) 
and ISewilo. et al. I (120041 ) have not yet been confirmed with high resolution ratio data. We 
exclude 11 UC sources, 2 compact sources and 2 diffuse sources that do not have significant 
continuum emission. 

We remove a. 
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(120061 ). Both of these recent catalogs compute the spectral index of SNR 
candidates using 20 cm and 90 cm VLA continuum data and rely on the ant i correl ation 
of SNRs with infrare d emission (8 nm MSX data in the case of iBrogan et al. I (120061 ) and 
21 yum MSX data for iHelfand et al. I (120061 )). We e xclude 13 SNRs from our sample tha t 



were found in these catalogs, as well as one found in iGaensler. Gotthelf. fc Vasisht I (119991 ). 
There are a comparable number of sources that are spatially coincident with SNRs, but 
which have a strong IR component. We believe these are Hll regions in locations that have 
produced multiple generations of stars. These sources are retained in our sample. We also 
remove an additional six sources that do not have infrared emission since they most likely 
are non-thermal. 

Finally, we determine if the classification (UC, compact, diffuse) and position of each 
source are correct, and remove duplicate sources. We require UC Hll regions to have small 
(< 1') bright knots of continuum emission, compact regions to be larger bright continuum 
sources, and diffuse regions to have faint extended continuum emission. There are many 
cases where an UC region was mistakenly identified as a compact or diffuse region in the 
L89 and L96 catalog s. This misident i ficatio n is due to the fact that UCs are unresolved 
with the 2^6 beam of lAltenhoff et al. I (119791 ) from which L89 and L96 drew their positions. 
We exclude 69 compact and diffuse nebulae whose positions are coincident with UC regions. 
Many of the UC H ll regions found in the GRS are in large complexes with individual UC 
components separated by angular distances less than the 46" GRS beam. We treat these 
complexes as single UC regions because they share a common molecular gas clump at the 
GRS resolution. 

Table [1] lists the 38 Hll regions we cull from our sample because of the criteria just 
described. Table [T] gives the source name, the reason for exclusion from our sample, and the 



^Available at 



http://www.mrao.cam.ac.Uk/surveys/snrs/l 
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reference if the source is a known SNR. Source names identify the nebular type: UC ("U"), 
compact ("C"), or diffuse ("D") . This source name convention will be followed throughout 
this paper. Figure [T] shows the longitude-velocity position of our nebular sample. The 
symbols in Figure [T] indicate the nebular type: UC (small filled circles), compact (medium 
filled circles), or diffuse (large open circles). Table [2] gives the properties of the nebulae in 
our sample. Listed are the source name, its position in Galactic and equatorial coordinates 
and its RRL velocit}!^! with its la error. Altogether our 301 Hll region sources probe 266 
unique directions since 33 nebulae have RRL emission at several different velocities that 
presumably originates from physically distinct nebulae located along the line of sight. For 
some Hll regions, we change the classification based on the morphology of the VGPS and 
MAGPIS radio continuum emission. We also change the position of a few H ll regions based 
on this continuum emission if the position is obviously incorrect. All such changes are noted 
with footnotes in Table [21 

4. Finding Molecular Gas Associated with Galactic H II Regions 

To characterize the properties of molecular gas associated with Galactic Hll regions 
we must establish reliable morphological correlations in (/, b, F)-space between the ^^CO 
gas and the nebulae. The GRS is a large (~ 10 Gbyte), dataset that contains very complex 
^^CO emission line structure. The (/, 6, Vlsr) structure of the molecular gas near the nebulae 
can be very complicated. Any given line of sight frequently contains multiple emission lines. 
There probably is no single algorithm that can uniquely give reliable CO/Hll regoim (/, b, 
V) morphological correlations. We therefore use a suite of software tools to analyze the GRS 
data cubes in a variety of ways. 

4.1. Analysis Software 

To maximize the power and flexibility of our analysis we wrote a large suite of IDL pro- 
cedures for spectral analysis rather than use any of the standard single-dish radio astronomy 
software packages. We used the single-dish radio astronomy TMBIDL0 software as our start- 
ing point. TMBIDL was originally written to analyze NRAO Green Bank Telescope data. 



•^The RRL velocities here are in the kinematic Local Standard of Rest (LSR) frame using the radio 
definition of the Doppler shift. The kinematic LSR is defined by a solar motion of 20.0 kms^^ toward {a, S) 
= (f8'\+30°)[f900.0]. 

"^Written by T. M. Bania and available at [http://www.bu.edu/iar/research/dapsdr/ 
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(TMBIDL was the inspiration that led to the NRAO GBTIDlI software.) The TMBIDL 
software emulates and improves upon many of the features of the NRAO UniPops analysis 
program. It includes Gaussian and polynomial line fitting, data visualization, data manipu- 
lation, etc. The TMBIDL code can easily be modified to analyze data from any single dish 
radio telescope. 

We wrote additional IDL procedures to interface and analyze GRS data within the 
TMBIDL environment. TMBIDL was created to analyze single spectra, so we added (l-V), 
(b-V), and (l-b) mapping tools to better visualize the GRS {l,b,V) data cubes. The basic 
visualization is a normalized contour map of the CO emission. Using these tools, we found 
that the morphology of the CO emission at the positions of the H ll regions was often very 
complex. 

To gain further control over the visualization of the GRS data, we wrote GUI-based 
softwarqjto analyze images extracted from the GRS {l,b, V) FITS data cubes. This software 
provides powerful GUI tools to extract, image, and analyze subcubes for each sample Hll 
region. For our analysis we imaged various quantities for an {l,b) zone surrounding each 
nebula. The user can, for example, easily modify the velocity channel whose ^^CO line 
intensity is being imaged over the mapped region. One can also quickly create and display 
integrated intensity CO maps, Wqo [K kms~^], as well as arbitrarily vary the velocity 
range of the integration, AV. (This is done with a kernel based algorithm that is fast and 
efficient.) Sub- images and regions can be created for any image, saved and then reloaded 
at any time. The GUI uses Dsd3 syntax to define regions. These regions can have all the 
basic DS9 shapes. We added a "threshold" region that selects all contiguous pixels above a 
user-defined threshold level that surround a given image pixel. Using this thresholding tool 
one can identify and analyze arbitrarily complex morphologies. 

This GUI-software is fully integrated with TMBIDL. The user can, for example, export 
the spectra within any region to TMBIDL for spectral analysis. The ability to scan quickly 
through velocity channels, create integrated intensity images, select pixels and fit Gaussians 
to spectra — all within a single application — is extremely powerful. 



^See http://gbtidl.sotirceforge.net/ 

^Available for download at http://people.bu.edu/andersld/ 
''Available for download at http://hea-www.harvard.edu/RD/ds9/ 
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4.2. Correlation of Molecular Gas and H II Regions 

Our goal is to find a morphological coincidence in {l,b,V) - space between the CO 
gas and the Hll region. After a coincidence is established, we want to characterize the 
molecular gas using spectral fits to the ^^CO emission. The UC positions are in general 
known to an accuracy greater than the 22" GRS pixel spacing. For the compact and diffuse 
Hll regions, our positions are accurate to a few arcminutes. The RRL LSR velocities are 
accurate to ~ O.lkms"^ (from Gaussian fits). Although the {l,b,V) position of each nebula 
is accurately known, establishing a robust set of criteria for identifying a real molecular/Hll 
physical association is nontrivial. 

The CO emission maps of Hll regions can be quite complex due to CMC structure 
and PDR/ionization front interactions. For example, one expects the molecular and ionized 
gas velocities to diverge as an Hll region evolves. Once an OB star forms within a CMC 
its ionization front (IF) expands rapidly at first, reaching the nebular Stromgren radius in 



5^, 10^ yr. The IF pushes the surrounding CMC molecular gas outwards. iDyson fc Williams 



(119971 ) show that as the IF expands, it rapidly slows until it is expanding at ~ lOkms"^ 
when it reaches the Stromgren radius. Since Hll regions are generally older than 10^ yr 
(with the possible exception of some UC Hll regions), the maximum difference between the 
RRL and the associated molecular gas should be < lOkms"^. 

Our nebular sample has H ll regions of different ages and thus should show evolutionary 
effects. Diffuse Hll regions should be older than the UC nebulae, and thus have had more 
time to evolve away from and displace their natal clouds. The molecular gas in diffuse 
H II regions should show a weaker association or may not be present at all. We expect the 
majority of UC an d compact Hll regions to be associated with a molecular clump (see, e.g.. 



Kim fc Koo II2OO3I ). Because of these complications, our CO/Hll analysis is comprised of a 



series of distinct investigations. 



4-2.1. Single Position Spectrum Analysis 



We first examine the GRS ^'^CO spectrum at the nominal position of each Hll region. 
Most previous studies of the molecular component o f H ll regions were made using single 
pointings, e.g.. IWhiteoak. Otrupcek. fc Rennie I (119821): Russiel fc Castets (2004). All report 
that the majority of H ll regions have associated CO. IWhiteoak. Otrupcek. fc Rennie I (119821 ) 
found in their survey of ^^CO emission from Southern H ll regions that molecular gas within 5 
kms~^ of the RRL velocity had large line intensities, and therefore was probably associated 
with the Hll region. In their analysis of Southern compact Hll regions using both ^^CO and 
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^^CO, iRusseil fc Castets I ( 2004J ) argued that 10 kms ^ is a better criterion for determining 
a molecular /H II association. 

To make a single pointing CO/Hll comparison we calculate an average ^'^CO spectrum 
at the position of each Hll region by convolving the GRS datacube (which is oversampled 
in angle) with the FCRAO telescope beam (HPBW = 46"). We then search this average 
spectrum for a ^^CO emission line peak at the Hll region RRL velocity. Specifically, we look 
for emission above 0.5 K brightness temperature, Tmb, and within ±2.5kms~^ of the Hll 
region LSR velocity. 

Our brightness temperature limit is chosen to be well above the GRS noise. The GRS 
data have a typical RMS sensitivity of cr(TMB) = 0.27 K. The beam convolved average 
spectrum has a factor of ~ 2 decrease in noise compared to a single GRS position. Thus 
our 0.5 K search criterion is a ~ 3cr limit. 

Using the 0.5 K and ±2.5kms~^ criteria, only 52 % of the nebular sample shows a 
CO/Hll association. Repeating this procedure with the same intensity requirement, but 
with velocity ranges of ±5.0kms~^ and ±7.5kms~^, we find that, respectively, 70% and 
79% of H II regions meet these criteria. Certainly increasing the velocity range further still 
will yield a greater number of nebulae matching the association criteria, but relaxing the 
association definition in this way also increases the chance of a misidentification and the 
possibihty of blending multiple velocity components. 



4-2.2. Wco Integrated Intensity Map Analysis 

We use the GRS data to make an (/, b) contour map of the ^'^CO integrated intensity, 
Wco (K kms~^), in order to provide information about the spatial distribution of the molec- 
ular gas. The main weakness of the single pointing method of searching for CO/Hll region 
associations is that sources with molecular gas offset from the nominal position of the Hll 
region are not counted as detections. For each nebula we use TMBIDL to make normalized 
Wco contour maps that are I x b = 12' x 12' (~ 30 x 30 GRS pixels) in size. The map 1^00 
is calculated by integrating each spectrum over 15kms~^ centered at the RRL velocity. The 
map peak Wco is used to normalize the Wco value for each pixel. The final product is a 
normalized contour map for each nebula in the sample. 

We then search each map for Wco peaks and note the distance from the 80% peak Wco 
contour to the nominal Hll region position. We define any source where this distance lies 
within the H ll region positional error bars to be a positive detection and any source where 
this distance is just outside of the error bars (roughly twice the positional uncertainty) to be 
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an ambiguous detection. All sources not meeting these positional criteria are deemed to be 
non-detections. We find that 70% of our sample nebulae show positive detections, 14% have 
ambiguous detections, and 16% show no correlation between the H ll region position and the 
^^CO emission. Somewhat suprisingly, adding the molecular spatial distribution information 
to the CO/Hll association criterion did not add significantly to the detection rate. 

4-2.3. CO (l,b,V) Data Cube Analysis Procedure 

Clearly, the molecular emission in the GRS is complicated and difficult to characterize. 
These experiments in establishing a CO/Hll region association demonstrate the need for a 
more sophisticated analysis. We use the §4.11 GUI software to search the GRS ^^CO data 
cubes in (/, b, V) parameter space. For each nebula we make a series of {I, b) Wco images and 
search for CO/Hll region associations. We follow the four step iterative procedure described 
below and illustrated in Figure [2J 

1. We first find the velocity range of the molecular emission associated with the Hll 
region. We examine single velocity channel images (position-position (/, b) maps) 
centered at the nominal position of the Hll region at its RRL velocity with overlaid 
VGPS 21cm continuum contours. We scan through these single channel images over 
ilOkms"^ of the source RRL velocity, searching for the channel where the molecular 
emission at the position of the H ll region has the highest intensity. If we are able to 
identify a molecular clump near the position of the H ll region, we extract the spectra 
from the voxels near this molecular clump. We fit a Gaussian to the (unweighted) 
average spectrum of this extracted emission and record the line center and FWHM 
of the emission line. Frequently, the molecular emission at the Hll region position is 
either absent or has a morphology that is difficult to characterize from single channel 
maps. For ~ 25% of our nebulae we are unable to make a molecular gas association 
from these single channel maps. 

2. Next we make an integrated intensity map, Wco-, by summing the intensities at a given 
{l,b) over the range of velocities found in step (1). If the source has an unambiguous 
CO/Hll association, we calculate Wqq over the velocity range V ± AV^/2, where V 
is the step (1) Gaussian line velocity and AV^ is its FWHM line width. If the source 
has an ambiguous step (1) association, we calculate Wco centered at the source RRL 
velocity over the range ±10kms~^. 

3. We then use the Wqq image created in step (2) to find pixels with molecular emission 
associated with the Hll region. We find the brightest emission near the Hll region 
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in the Wqo image and select all contiguous pixels that have values above a threshold 
determined independently for each Hll region. The threshold is varied until a small 
number of pixels are selected — typically 20 to 30. The exact number is set by the 
molecular clump's intensity profile. Small clumps with a sharply peaked intensity 
profile have fewer pixels, whereas larger clumps with a "plateau" of emission have 
more pixels. This method seeks to isolate single clumps in order to preserve the line 
peak intensity and minimize blending of multiple velocity components. 

4. All the GRS spectra within the (/, b) region selected in step (3) are used to calculate 
an average ^^CO spectrum. We then fit the fewest possible number of Gaussian com- 
ponents to the spectrum in order to maximize the intensity of the peak CO emission. 
The majority (~ 80%) of our sources are adequately fit with only one Gaussian com- 
ponent. Many nebulae, however, do not have clean Gaussians profiles, but rather show 
structure that often suggests a fainter, wider hne superposed on a brighter, narrower 
line. 

Our analysis procedure is summarized in Figure [2] for the Hll region U43.24— 0.05. 
Panels A and B depict step 1. Panel A shows a single channel GRS map at the velocity 
nearest the RRL velocity where the molecular clump has the highest intensity. The cross 
marks the nominal Hll region position. For clarity we have not shown the VGPS continuum 
emission contours. The black box shows the (/, b) positions of the voxels from which we 
extract spectra to produce the average spectrum shown in Panel B. The vertical line in 
Panels B and D marks the RRL velocity. Using the velocity range of the emission line shown 
as dashed lines in Panel B, we create the integrated intensity image shown in Panel C (step 
2). This image has been smoothed with a 3 x 3 Gaussian kernel. The black outline in 
Panel C shows contiguous integrated intensity values above a threshold (step 3). This region 
and its threshhold is defined by visual inspection of the image. This is the ^^CO emission 
we deem to be associated with the Hll region. Finally, we extract the spectra from these 
(/, b) positions to produce the average spectrum shown in Panel D (step 4). For this source, 
our method separates the two emission lines that are blended in the panel B spectrum and 
preserves the peak line intensity. 

In principle, we could iterate our analysis to locate and fit the molecular emission more 
accurately. Using the Panel D Gaussian fit, we could create a new integrated intensity image, 
define a new CO/Hll association region, and fit a Gaussian to this new average spectrum. 
We did this for 10 test cases and found only only minimal changes that were not significantly 
different from a single pass analysis. 

This procedure has many advantages. We are able to characterize the spectral prop- 
erties of molecular gas distributions that have arbitrary morphologies. We minimize our 
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assumptions at every stage of the process. By first examining the {l,b) images of each source 
at individual velocity channels, we limit false detections that may arise from integrated in- 
tensity images containing multiple velocity components blended together. By extracting the 
spectra from regions defined in the integrated intensity images, we make no assumptions 
based on the visual appearance of the molecular emission at individual velocity channels. 
By using a variable threshold to select pixels, we are able to characterize molecular struc- 
tures with arbitrary morphology. This threshold definition ensures clean spectral fitting; 
the spectra are not contaminated by adjacent pixels that would lower the line intensity and 
might increase the width of the fitted Gaussian line. Finally, since we analyze an average of 
many spectra, the lines we fit have a much greater signal to noise ratio than a single pointing 
spectrum. The Gaussian fit uncertainties in the spectral line parameters we derive are thus 
minimized. 

After we establish a CO /H ll region association, we characterize the angular size of the 
molecular cloud by fitting an ellipse to the Wco spatial distribution. Since the association 
defined in Panel C of Figure [2] is a Wco threshold value that is determined independently 
for each source, the size of an ellipse fitted to this zone would have little physical meaning. 
We therefore define a new threshold selected region that is uniform for the entire nebular 
sample. This new region is defined by a threshold set to 80% of the Wco peak inside the 
original association zone. The fitting algorithm calculates the semimajor and semiminor 
axes using the "mass density" of pixel locations: the clustering of pixel locations along the 
/ and b directions. The fitted ellipse to this region is thus a uniform estimate of the size of 
the molecular gas associated with the H ll region. 

The properties we derive for the CO/Hll region associations are summarized in Table 
131 For each nebula we list the parameters of the ellipse and Gaussian line fits. All errors are 
la uncertainties. Given are the ellipse centroid in Galactic coordinates, its size, semi-major, 
a, and semi-minor, b, axes, together with the position angle measured from north toward 
increasing Galactic longitude. The ellipse size is the geometric mean diameter, 2^/ab. If there 
are multiple spectral components in the Gaussian fit only the properties of the brightest are 
listed. We also use only this brightest component in our subsequent analysis. The fitted line 
parameters given are the center velocity, intensity (in main beam brightness temperature 
units), and FWHM line width, AV. Table [3] also lists some additional nebular properties 
derived below in ^ the CO excitation temperature, the ^^CO column density, and the 
nebular confidence parameter, CP. 



-14- 



5. Discussion 

Our analysis here provides a large sample of molecular cloud/Hll region associations 
wliose physical properties are well characterized. In ^we show that the molecular emission 
is often morphologically complex and offset from the nominal H ll region position and RRL 
velocity. We find that the traditional single pointing analysis does not reliably detect the 
molecular components of Hll regions. It is very difficult to make a CO/Hll association with 
any confidence without an image produced from a (/, b, V) data cube that spans the entire 
Hll region/PDR/molecular cloud interaction region. 

To distinguish sources with an unambiguous molecular gas component at the H ll region 
position and velocity from those with less robust molecular gas associations, we assign a 
confidence parameter, CP, to each source ranging from A to E. Our qualitative criteria 
for the confidence parameter are as follows: A: no ambiguity in position or velocity, the 
molecular gas coincident with radio continuum or shows clear signs of interaction; B: either 
offset somewhat in position or velocity, or in a complex region of molecular gas; C: either 
offset in position or velocity or in a complex region, but fainter than a B source and offset 
further; D: diffuse emission near the correct position and velocity, but uncharacterizable 
due to low intensity or ambiguous morphology; E: nothing at all apparent in position and 
velocity. 

Figures [3], H] and [5] give representative examples of our confidence parameter classifi- 
cation. These 20' x 20' ^^CO integrated intensity images span the range of molecular gas 
morphologies surrounding our sample of H II regions. The images are grouped by their confi- 
dence parameter classification: Figure [3] shows CPA sources. Figure H] shows CP B sources 
and Figure shows CP C, D, and E sources. Because the GRS is oversampled (22" pixels at 
46" HPBW), we increase the signal to noise by smoothing the images with a 3 x 3 Gaussian 
filter. 

At the top of each image we list the H ll region name, the ^^CO line velocity ( km s~^) and 
FWHM line width (kms^^), AV, and the source's confidence parameter. These Gaussian 
fitted line parameters are used to calculate the source's ^^CO integrated intensity, Wco(K 
kms~^), for the velocity range V ± AV/2. At the bottom of each image we give the Wco 
calculated from the fitted line: Wco = 1.06 Tmr Af , which is accurate for Gaussian line 



shapes. Using the procedure described in ISimon et al. I (1200 ll ). we also give an estimate of 
the H2 column density, A^(H2)[cm~^]. 

The grey scale image shows the Wqo distribution. Normalized contours are drawn at 
83%, 67%, 50%, 33% and 16% of this maximum. (One can get quantitive Wqo values using 
the scale bar at the right.) For Hll regions where we could not associate molecular gas, the 
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contour lines are dashed rather than sohd and, of course, no hne parameters are given. 

A bold cross marks the nominal position of the Hll region. Plotted in bold is the 
fitted ellipse described in §4.2.31 If there are other H ll regions in the field, they are marked 
with thinner crosses. The cross arm lengths correspond to the beam size used to make the 
measurement of the Hll region (3' for L89; 9' for L96), except for the UC regions where the 
cross arm lengths are set to 1'. The GRS beam (HPBW = 46") is shown in the lower left 
corner of each image. Shown in grey in these images are VGPS 21cm continuum contours. 
Tickmarks on these contours point downhill, towards decreasing 21cm emission. 

We are able to establish a highly confident {CP A and B) CO/Hll region association 
for 62% of the nebulae in our sample. Relaxing the confidence criterion to CP A, B and C 
sources gives CO/Hll associations for 84% of our nebular sample. Histograms of the number 
distribution of confidence parameter values are shown in Figure E] for UC (open), compact 
(hatched), and diffuse (gray) nebulae. The top left panel is the stacked histogram of the 
distribution; the top line represents the entire sample. For example, there are a total of 113 
CP A sources: 78 UC, 29 compact, and 6 diffuse Hll regions. Clockwise from here are the 
individual histograms for the UC, compact, and diffuse nebulae. (Most subsequent figures 
will follow this display format.) As is clear from Figure O the UC sample has the greatest 
number of high confidence CO/Hll region associations. 

Some Hll regions appear to have no molecular gas associated with them. Of our 301 
sources, 14 (5%) are classified as E sources; these nebulae show no ^'^CO emission whatsoever. 
Thirty four sources (11%) are classified as D and therefore have only diffuse emission at the 
correct position and velocity. Five of our E sources and 15 of our D sources have multiple 
RRL velocities along the same line of sight. In these cases, one RRL velocity is probably 
from the Hll region of interest while the other is likely from a nearby Hll region. Eleven 
of the 12 UC regions with confidence parameter values of D or E have multiple velocity 
components. 

Two UC H II regions are worth mentioning individually because of the nature of their 
molecular associations. The UC Hll region U33. 13+0.09 has a very large velocity offset 
between the molecular m a terial and the RRL velocity. L89 find a RRL velocity of 93.8 



kms ^ and lArava et al. I (120021 ) find a RRL velocity of 87.4 kms \ Our ^^CO velocity 



of 75 kms ^ is in agreement with the CS velocity found by iBronfman. Nyman. fc May 



(119961 ). The morphology and linewidth of the molecular emission suggest that the molecular 
gas is associated with the Hll region. The UC Hll region U21.42— 0.54 also has a compact 
molecular clump at the correct (/, b) position. The ^^CO velocity of this clump, 54 kms~^, 
is 16 kms~^ offset from its RRL ve locity of 70 kms~^. This source was not detected 



by [Bronfman. Nyman. fc May I (119961 ). We assign a confidence parameter value of C to 
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U21.42— 0.54 because of the extreme velocity offset. 

The lack of CO/Hll region associations in ~ 20% of our sample is not entirely unex- 



Blitz. Fich. & Stark 


( 


1982) 


studied. 


Russeil & Castets 



( I2OO4J ) f ound no association with CO in ~ , 2 0% o f their sample of Southern compact Hll 



regions, 
monia in ~ 



Churchwell. Walmsley. fc Cesaroni I (jl990l ) do not detect the dense gas tracer am- 



30% of a sample of 84 UC H ll regions. 



Our large sample of Hll regions contains nebulae spanning a range of evolutionary 
stages. The UC regions, being young, are more likely to lie within their natal molecular 
clouds where the density of ^^CO should be high. As the Hll region evolves, it will dissipate 
the gas. We should see evidence for this in the compact and diffuse regions. We therefore 
expect UC regions to be associated with small sizes, high excitation temperatures and column 
densities, large line widths, and bright line intensities. Diffuse regions should have large sizes 
and smaller values for the other quantities compared to the UC regions. Compact regions 
should lie in between these two extremes. 

The lack of molecular gas in an H ll region can certainly be an evolutionary effect. These 
nebulae may represent an older population of Hll regions that have had time to become 
displaced from the gas in which they formed due to a variety of mechanisms including, stellar 
winds, ionization fronts, high stellar space velocities, etc. One observational consequence 
of this scer iario might be a bubble morphology in the molecular gas, seen as a ring in 
projecti on. Churchwell et al. I ( 2006 ) found a large population of bubbles in the GLIMPSE 
survey (IBenjamin et al. II2003I ). Observed at mid-infrared wavelengths, GLIMPSE boasts 
an angular resolution 10 times that of the GRS and is therefore a better diagnostic tool for 
locating bubbl e features. We indeed see bubbles in ^^CO in most H ll regions in the GRS for 

(I2OO6I ) also identified bubbles. These regions are usually 



sources where Churchwell et al. 



classified as D or E sources because the gas has been pushed far away from the center of the 
continuum emission. Further investigation of this topic will be the subject of a future paper. 



5.1. Properties of Molecular Cloud/H II Region Sources 

Here we focus on the 253 nebulae with positive CO/Hll associations, the subset of our 
Hll region sample with confidence parameter values of A, B or C. The (/, V) distribution of 
these nebulae is shown in Figure [7] for CPA (large filled circles), B (small filled circles), and 
C (small open circles) sources. Fully 90% of UC and compact nebulae have CP A, B, or C 
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quality CO/Hll associations whereas only 64% of the diffuse nebulae do. Diffuse Hll regions 
are probably older on average than either UC or compact nebulae so their significantly lower 
associate rate with molecular gas provides the first hint of evolutionary effects in our sample. 

Table H] summarizes the mean properties we derive here for this sample. Listed are 
the mean and standard deviation (la) for each quantity. This information is given for the 
entire sample and also for various subsets of it: UC, compact, and diffuse Hll regions as 
well as sources of CP A, B and C. Table H] lists the number of Hll regions in each category, 
the absolute value of the velocity difference between the ^^CO molecular clump and the 
RRL Vlsr, the line intensity (main beam brightness temperature), the FWHM line width, 
the size of the associated molecular clump, the CO excitation temperature and the ^^CO 
column density. 



We estimate the average CO column density towards our sources using the lRohlfs &: Wilson 



(119961 ) analysis: 



N{^^GO)[cm-''] = 2.42 x 10 



14 



1 -exp(-5.29)/Te, 



where Tex is the excitation temperature and ris is the optical depth of the ^^CO line. We 
assume the ^^CO emission is optically thin and use the Gaussian fit ^^CO line parameters 
to find the optical depth integral in EqO Both the optical depth and Equation [1], however, 
depend on the excitation temperature, T^^. 



We use the UMSB ^^CO J 



survey (ISanders et al. Ill986l ) to estimate Tex for 



each source. We assume ^^CO is optically thick and use the radiative transfer equation for 
the J = 1 — > transition to calculate the excitation temperature from the observed main 
beam brightness temperature, Tjjjg, of the ^^CO hne: 



Tex = 5.5 /ln( 1 + — ^^^^ ). (2) 

Eq. [2]holds as long as: (1) the ^^CO and ^^CO emitting gas is in LTE at the same excitation 
temperature; (2) this gas fills the same volume without clumping; and (3) there are no 
background continuum sources. (For our nebulae the Hll region continuum is subtracted 
when the spectral baselines are removed.) 

For each nebula we first compute an average spectrum from the UMSB survey datacube 
in exactly the same way as we did for the GRS data. We use the same {I, h) positions 
from the identical threshold selected region for this average. We then fit Gaussians to these 
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average spectra. If we fit multiple Gaussians to the GRS data, we attempt to fit the same 
components to the UMSB data. The spectral resolution of the UMSB survey is lkms~^, so 
this was not always possible as lines resolved in the GRS are blended in the UMSB survey. 

Because the pointings of the UMSB survey are further apart than in the GRS (180" com- 
pared to 22"), the UMSB survey underestimates the ^^CO emission for the small molecular 
clumps found in the GRS. For a given velocity, each pixel in the UMSB survey represents 
~ 70 pixels in the GRS. Each threshold selected region contains ~ 20 GRS pixels on average, 
so for the majority of our sources we use only one UMSB pointing to estimate the excitation 
temperature appropriate for the ^^CO emission. Furthermore, this pointing can be as far 
away as ~ 120" from the GRS position. 

The distribution of excitation temperatures we derive is shown in Figure [3 All H ll 
regions in our sample have very similar excitation temperatures near the standard 10 K 
value assumed for molecular clouds. We expected the excitation temperature of the UC 
regions in particular to be higher than this standard value as the CO gas is nearer to the 
exciting star. That we do not see hotter temperatures associated with younger regions is 
probably due to the effect of the undersampling of the ^^CO emission by the UMSB survey. 
The mean ^^CO to ^^CO ratio for molecular clumps smaller than the beam spacing of the 
UMSB survey, 3', is ~ 2 while this ratio is ~ 3 for molecular clumps larger than 3'. Our 
excitation temperatures, and hence column densities, are therefore lower limits. Because of 
their small size, UCs are affected more by the difference in sampling between the GRS and 
the UMSB survey. The UCs also suffer from beam dilution which will lower the inferred 
excitation temperature. 

Nevertheless, we use these excitation temperatures to compute the ^^CO column density 
for each source using Eq. [H These column densities are listed in Table [2J We then estimate 
the H2 column density. 



7V(H2) 



•12 



CO" 

13C0 



12C0 



(3) 



by assuming constant values for these abundance ratios. Following ISimon et al. I (120011 ) we 
adopt a i^CO/i^CO ratio of 45 and a Hs/^^CO ratio of 8 x 10"^ Our ^^CO and H2 column 
density results are summarized in Figure [91 As expected, the UC nebulae have on average 
the highest column densities and the diffuse nebulae the lowest. Too, the CPA sources have 
on average over twice the column density of the CP C nebulae. 



We find that the CO gas has on average only a small velocity offset from the Hll region 
RRL velocity. Figure [10] shows the difference between the velocity of the CO gas and the 
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RRL velocity. There is no difference in velocity offset between the various types of Hll 
regions. The Gaussian fit to the entire distribution is centered at 0.4 kms'^ with a FWHM 
of 8.5 kms~^, whereas the mean of the distribution is 0.2 ± 3.8 kms^^. The fact that the 
distribution is centered at zero velocity offset is to be expected for a RRL selected sample 
of Hll regions and an optically thin tracer such as ^^CO. 

This result is in contrast to optically selec ted samples where a positive velocity offset 



was found (e.g., iFich. Dahl. fc Treffers I Il990l ). Optical samples choose specific CO/Hll 



region line of sight geometries. Face-on or edge-on blister sources such as the Orion nebula 
and M17 should dominate these samples. In our sample there is no preferred radial location 
for Hll regions within molecular clouds. Because we do not know the CO/Hll geometry 
with respect to the line of sight, the absolute value of the CO/RRL velocity difference will 
be a direct measure of any systematic velocity offset between the molecular and ionized gas. 
Table m therefore lists the mean absolute value of the CO/RRL velocity difference; it shows 
that there is a ~ 3kms~^ average flow velocity between the molecular and ionized gas for 
the nebulae in our sample. This value is independent of the type of H ll region, but increases 
slightly as the CP value decreases from A to C. 

Figure [11] shows the distribution of ^'^CO line intensities for our sample. Based on 
the evolutionary model of Hll regions the natal cloud is gradually dissipated by photo- 
dissociation, photo-ionization, and expanding motions. We expect the ^^CO density to 
decrease as the region progresses from UC to compact and then, finally to diffuse. Assuming 
^^CO is optically thin, or at least marginally so, the higher densities found in UC regions 
would lead to higher line intensities, whereas compact sources should show lower line in- 
tensities, and diffuse sources the lowest. This hypothesis is only partially borne out as UC 
and compact regions share the same distribution, averaging 5.21 ± 0.23 K and 4.96 ± 0.24 K, 
respectively. Diffuse regions do show lower line intensities of 3.32 ± 0.19. The errors quoted 
here are the standard errors of the mean, s.e.m. = a/^/N. 

Figure [12] shows the distribution of ^'^CO line widths for our sample. The Gaussian fit 
to this distribution is centered at 4.0 kms~^ with a FWHM of 3.2 kms~^. The mean of this 
distribution is 4.20 ±0.09 kms~^. We expected UC Hll regions to have significantly broader 
lines than compact Hll regions because the molecular gas is closer to the exciting star and 
the outflows should be stronger. Once again, this is not the case: the UC and compact 
distributions are very similar, averaging 4.43 ± 0.13 kms"^ (s.e.m) and 4.23 ± 0.14kms"^ 
(s.e.m.), respectively. UC and compact regions do have broader lines than the diffuse regions 
which average 3.56±0.20 km s~^ (s.e.m.) . This suggests that the central star(s) may no longer 
be significantly heating molecular gas near diffuse H ll regions. 



This distribution of line widths is comparable to that found by 



Russeil &: Castets 



(120041 ) 
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in their single pointing survey of southern H ll regions. They find that the ^^CO J = 1 ^ 
fine has an average fine w idth of 3.7 kms~^ w ith a standard deviation of 1.9 kms^^. In a 



study of UC Hll regions, iKim &: Koo I ( 20031 ) find an average fine width of 6.8 kms ^ in 



the ^^CO J = 1 transition. Their calculation of line width, however, was based on the 
average spectrum over the entire map area, w hich was as large as 30' x 40'. For the 8 UC 



regions in our sample that iKim fc Koo I ( 20031 ) study, we measure a line width of 5.0 kms ^ 



whereas they measure 7.2 kms ^. Using the same large areas to compute the line widths 
for these sources, we find an average line width of 7.3 kms~^. 

The angular size distribution of the H llCO sources is shown in Figure [131 These sizes 
are defined as the geometric mean of the major (2 x a) and minor (2 x b) axes of the fitted 
ellipse, 2\/ab. (See §4.2.11 for our ellipse fitting procedure.) There are 6 sources that have 
sizes greater than 10' that are not plotted in Figure [131 These sources are invariably clumps 
in a large region of extended molecular emission, which makes our method of determining 
the angular size unreliable. 

The molecular clumps associated with UC regions do show the smallest sizes, as ex- 
pected, averaging 1^7 ± Ofl (s.e.m.). Compact Hll regions are slightly larger, averaging 
2'2 ± O'l (s.e.m.). The average size of the molecular gas associated with diffuse Hll regions 
lies in between that of UC and compact Hll regions, averaging 1'9 ± O'l (s.e.m.). We have 
removed the 6 sources with sizes greater than 10' from the statistical analysis. The molec- 
ular gas around many diffuse regions is fragmented, which leads to the small angular sizes 
we measure. Since we only associate a single molecular clump of contiguous pixels with 
each nebula, for diffuse Hll regions we probably have not characterized all the associated 
molecular gas. 



5.2. Comparison with GRS Molecular Clumps 

The properties of molecular clumps in the GRS were analyzed down to size scales of 
~ 1' (Rathbor ne et al. 2008, in preparation). The contiguous pixel finding algorithm 
CL UMPFIND dwiUiams. de Geus. fc Bhtz Ill994[ l was used to locate Giant Molecular Clouds 



(GMCs) within the GRS. Then, by altering the size threshold in CLUMPFIND, the clumps 
within the GMCs were identified and characterized. The distribution of peak intensities and 
line widths for these clumps shows a Gaussian core with an exponential tail at high values of 
each parameter. The break points where the distributions turn over from being dominated 
by the Gaussian core to being dominated by the exponential tail are roughly 4 K and 2 
kms~^. By number the vast majority of these GRS molecular clumps have line intensities 
below 4 K and line widths below 2 kms~^. 
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The molecular gas associated with Hll regions has on average a greater line intensity 
and larger line width compared to molecular clumps in the GRS. We plot in Figure [H] the 
line intensity verses the FWHM line width for the molecular clumps associated with our Hll 
regions: UC (filled triangles), Compact (filled circles), and Diffuse (open squares) nebulae. 
The solid lines divide the plot into quadrants according t o the b reak points of the GRS 
clumps. A similar plot was used by IClemens fc Barvainid (119881 ) to show that the small 



clouds in their optically selected molecular cloud sample were cool and quiescent. 

The lower left quadrant in Figure [TH should be populated by cold quiescent clouds. These 
objects are neither making stars nor being externally heated. The upper left quadrant should 
contain a population of clumps that are heated externally. These clouds are warm (or have 
high column densities), but do not have the non-thermal motions that would be present 
if they posessed a central star. The upper right quadrant should contain molecular gas 
associated with embedded massive stars. The lower right quadrant should contain embedded 
protostars. These large line width objects are likely active sites of star formation, or near 
an active site. Thus the lower, < 4 K, part of the plot has a pre-stellar population whereas 
the upper part contains clouds that are affected by local massive stars. 

The UC and compact nebulae occupy the same region of Figure [TH they have similar 
molecular properties. The diffuse regions, however, have lower line intensities and line widths; 
they are similar to the general population of molecular clumps. These clumps are no longer 
being heated significantly by the ionizing star. 

Most of our sources with associated ^'^CO, 54%, lie in the upper right quadrant of Figure 
where large line intensities and broad line widths suggest active sta r formation. The vast 



majority of GRS clumps, as well as most of the molecular clouds in IClemens fc Barvainis 



(lOSSj), reside in the lower left quadrant. The bulk of our remaining nebulae, 42%, lie 



in the lower right quadrant. The UC, compact, and diffuse regions all have a significant 
population in this quadrant. It is tempting to think that these low line intensities are due 
to the decreased optical d epth of ^^CO but Figure [H] looks very similar to the same plot in 
J2OO4I ) made using the optically thick ^^CO J = 2 ^ 1 transition. 



Russeil &: Castets 



5.3. Are Molecular Cloud/H II Region Associations Real? 

Are these CO/Hll region associations really sources that are having a direct physical 
interaction between the Hll region and ambient molecular gas? The associations are based 
on morphological matches in position and velocity between the ^^CO gas, RRL velocity, and 
radio continuum emission. But correlation does not imply causality: these matches could in 
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principle be a coincidental juxtaposition projected on the sky into the same sohd angle by 
molecular clouds and H ll regions located at entirely different places along the line of sight. 

That we require a morphological match in (/, 6, V) - space places a severe constraint on 
a false positive association. Mere (/, h) coincidence is not enough; the velocity also needs to 
match. The kinematic distance ambiguity in the Inner Galaxy makes it possible for the H ll 
region and CO cloud to be at different line of sight positions despite having nearly identical 
radial velocities. But these are special places because only they share the same LSR velocity. 
Assessing the quantitative probabihty of a false positive association is beyond the scope of 
this paper. To our knowledge no one has yet done the detailed modelling this would require. 
One needs to know the Galactic distribution of the clouds which posits a detailed knowledge 
of Galactic structure. For a false positive association we require that there not be a cloud 
at the Hll region position, but that there be a cloud at the other kinematic distance. One 
thus needs to evaluate separately for each H ll region the line of sight distance derivative of 
the LSR velocity, dV/dr in order to assess the path lengths at the near and far kinematic 
distances that must be populated. This requires a detailed knowledge of Galactic kinematics, 
including streaming motions caused by spiral arms. With this information one might be able 
to estimate the probability of a false positive association. We probably do not know enough 
about either Galactic structure or Galactic kinematics to do this. 

The fact, however, that ~ 20% of the Hll regions do not have associated CO gas 
( §5.ip is evidence that suggests chance line of sight superpositions in (/, 6, V) - space are 
rare. Furthermore, Figure [TH provides strong support for the physical reality of our CO/Hll 
region associations. Our ^^CO clouds are not only near to the Hll regions in {l,h,V) - 
space, but their spectra also have the trademarks of star formation: bright lines and large 
line widths. This is in marked contrast to the spectral line properties of the vast majority 
of GRS clouds (see §5.2p . Only 1/3 of the GRS clumps have peak intensities > 4K whereas 
55% of our CO/Hll associated clouds do. Only 1/6 of the GRS clouds have line widths 
> 2kms~^; nearly all of our clouds, 96%, exceed this value. We conclude that most of 
the CO/Hll region associations must be nebulae with real physical interactions between the 
molecular and ionized gas. 

6. The GRS HII Region Catalog 

Our analysis here produced a catalog of Figure [3] type images and physical properties 
for the sample of 301 Galactic Hll regions. We created a websit^ to give everyone access 



http://www.bu.edu/iar/hii_regions 
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to this information. In addition to the images of nebular Wqo, this website has the average 
^^CO spectrum of each source as well as all the information found in Tables [2] and [31 We 
expect this website to be an evolving database compiling additional information about these 
nebulae as it becomes available. 



Future Work 



Our sample of CO /H ll region associations will enable many further studies of the prop- 
erties of star forming regions at all stages of their evolution. The most important parameter 
that is missing here is the distance to each nebula. Knowing the distance would enable us to 
derive the intrinsic physical properties of each nebula, establishing their physical sizes and 
turning column densities and line intensities into masses and luminosities. 



Anderson fc Bania I (120081 ) [AB hereafter] use H I absorption studies to derive kinematic 



distances toward all the Hll regions with associated molecular gas. All our nebulae are in 
the first Galactic quadrant, so their distances are degenerate due to the kinematic distance 
ambiguity. Using the fact that H I absorbs thermal continuum from the H ll region, AB use 



theVGPS 21cm Hi emission line maps to remove this degeneracy (see iKuchar fc Bania 
19941 ). This is a proven technique as there is sufficient residual col d H I associated with 



almost all GRS mo l ecular clouds to p roduce significant absorption (jjackson et al. 1 12002 



Kolpak et al. II2003I : iFlynn et al. II2004I ). We shall then use these distances to analyze this 
nebular sample and derive the physical properties o f the dust and gas (io nized, atomic & 
molecular). The completion of the Spitzer GLIMPSE ( iBenjamin et al. 1120031) and MIPSGA L 
(Carey et al. 2008 in preperat ion) surv eys, together with th e GRS (I Jackson et al. II2006I ). 
MAGPIS jHelfand et al. ibooeh . NVSS JCondon et al. Ill998h and the VGPS surveys enable 
for the first time a multi-wavelength analysis of the physical properties and evolutionary 
state of a large sample of inner Galaxy Hll regions. Due to our large sample size, we will be 
able to find examples of H ll regions at all evolutionary stages. 



8. Summary 

We analyzed the GRS ^^CO molecular gas associated with all known H ll regions covered 
by the GRS using multiple analysis techniques. Our sample includes 301 regions: 123 UC, 
105 compact and 73 diffuse H ll regions. We found that 80% of our H ll regions showed pos- 
itive molecular associations, with UCs having the highest association percentage and diffuse 
regions the lowest. About 5% of our sample showed no molecular emission whatsoever. We 
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hypothesize that some of these non-detections represent an older population of H ll regions 
where the molecular gas has been displaced from the central star or stars. We found that 
the molecular properties of UC and compact H ll regions are quite similar, with line widths 
averaging ~ 4kms~^ and ^^CO column densities of about 3.5 x 10^^ cm^^. The molecular 
gas associated with diffuse regions has properties more consistent with quiescent clouds. The 
molecular gas properties of our sample nebulae are consistent with an evolutionary sequence 
wherein small, dense molecular gas clumps associated with UC H ll regions grow into older 
compact nebulae and finally fragment and dissipate into large, diffuse nebulae. 

This publication makes use of molecular line data from the Boston University-FCRAO 
Galactic Ring Survey (GRS). The GRS is a joint project of Boston University and Five Col- 
lege Radio Astronomy Observatory, funded by the National Science Foundation under grants 
AST-9800334, AST-0098562, & AST-0100793. The National Radio Astronomy Observatory 
is a facility of the National Science Foundation operated under cooperative agreement by 
Associated Universities, Inc. 
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Table 1. Faux and Anomalous Sources 



Source 


Notes 


Reference 


C14.32+0.13 


SNR 


a 


D15.45+0.19 


SNR 


a 


D15.52-0.14 


SNR 


b 


U16. 58-0.05 


No continuum peak 




D17.23+0.39 


Probably an evolved star 


c, d 


U17.64+0.15 


No continuum peak 




C18. 64-0.29 


SNR 


a, e 


U19. 12-0.34 


No continuum peak 




U19. 36-0.02 


No continuum peak 




C19. 88-0.53 


No continuum peak 




C20.26-0.89 


No IR 




C20.48+0.17 


SNR 


a, e 


D21.56-0.il 


SNR 


a, e 


D22.04+0.05 


No continuum peak 




D2216-0.16 


Star 


c 


D22.40-0.37 


SNR 


e 


C22. 94-0.07 


No IR 




C23. 07-0.37 


No IR - Part of SNR? 




C23. 07-0.25 


No IR - Part of SNR? 




D23.16+0.02 


No continuum peak 




U23. 24-0.24 


No continuum peak 




D26.47+0.02 


WR or LBV 


f 


U26. 51+0.28 


No continuum peak 




C27.13-0.00 


SNR 


e 


C29. 09-0.71 


SNR 


e 


D29.55+0.11 


SNR 


g 


U30 42+0 46 


No continuum pG3,k 




D30.69-0.63 


No IR 




U30.82+0.27 


No continuum peak 




C30.85+0.13 


SNR 


e 


C31. 05+0.48 


SNR 


e 


D31. 61+0.33 


SNR 


e 


D31. 82-0.12 


SNR 


e 


U33. 24+0.01 


No continuum peak 




C45. 48+0.18 


No continuum peak 




U49. 67-0.45 


No continuum peak 




C50.23+0.33 


No IR 




U53. 63+0.02 


No continuum peak 




HBroean et al. 


1 I2OO6I) 




iBroean et al. 


1 l|2006l") Qow confidence detection^ 


ISteohenson 1 \l99i) 




tewok, Volk. i 


k Bidelman 1 lll997l') 




'iHelfand et al. 


1 J2OO6I) 





IClark et al~1 ll2003h 
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Table 2. Hll Region Source Sample 



Source 


1 


b 


RA(2000.0) 


DEC(2000.0) 


Vlsr 


Reference 


Comments 




n 


n 


(h m s) 


(O , „) 


(kms-i) 






D15.00+0.05a 


15.00 


+0.05 


18 17 41 


-15 52 50 


26.5 ± 1.6 


L96 




D15.00+0.05b 










63.5 ± 1.7 


L96 




D15.64-0.24 


15.64 


-0.24 


18 19 60 


-15 27 10 


61.8 ± 1.3 


L96 




C16.31-0.16 


16.31 


-0.16 


18 21 01 


-14 49 30 


49.5 ± 0.7 


L89 




C16.43-0.20 


16.43 


-0.20 


18 21 24 


-14 44 10 


44.5 ± 0.9 


L89 




D16.61-0.32 


16.61 


-0.32 


18 22 11 


-14 38 10 


44.9 ± 0.6 


L89 


a 


D16. 89+0.13 


16.89 


+0.13 


18 21 05 


-14 10 30 


42.3 ± 1.6 


L96 




D17.25-0.20a 


17.25 


-0.20 


18 22 59 


-14 00 50 


49.9 ± 1.4 


L96 




D17.25-0.20b 










96.5 ± 1.9 


L96 




U18.15-0.28 


18.15 


-0.28 


18 25 01 


-13 15 20 


53.9 ± 0.4 


L89 





''Changed nebular classification from compact 

Note. — Table [2]is published in its entirety in the electronic edition of the Astrophysical Journal Supplement 
Series. A portion is shown here for guidance regarding its form and content. 



Table 3. Properties of Molecular Cloud/HllRegion Sources 



Fitted Ellipse Parameters Fitted Gaussian Parameters 



Source 


1 


b 


Size 


Maj. X 


Min. 


PA 


V 


Tmb 




AV 






Af("CO) 


CP 




(°) 


(°) 


{') 


('X 


') 


n 


( kms^-'-) 


(K) 




( km s" 




K 


(xlOi'^cm-2) 




D15.00+0.05a 


14.93 


+0.02 


2.1 


1.6 X 


0.7 


+88.6 


25.85 ± 0.02 


6.03 ± 


0.02 


3.77 ± 


0.05 


12.7 


3.2 


B 


D15.00+0.05b 




























E 


D15. 64-0.24 


15.66 


-0.21 


2.0 


2.1 X 


0.5 


+70.6 


56.96 ± 0.06 


1.88 ± 


0.06 


1.86 ± 


0.09 


8.6 


0.4 


B 


C16.31-0.16 


16.36 


-0.21 


1.7 


0.9 X 


0.8 


+45.0 


47.56 ± 0.11 


4.17 ± 


0.11 


4.68 ± 


0.12 


11.0 


2.6 


B 


C16.43-0.20 


16.36 


-0.21 


1.7 


0.9 X 


0.8 


+14.9 


48.83 ± 0.02 


9.47 ± 


0.02 


2.74 ± 


0.04 


13.9 


3.9 


B 


D16. 61-0.32 


16.56 


-0.34 


4.8 


2.9 X 


2.0 


+26.3 


43.23 ± 0.01 


5.39 ± 


0.01 


4.44 ± 


0.04 


11.1 


3.2 


C 


D16.89+0.13 




























D 


D17.25-0.20a 


17.23 


-0.24 


3.9 


2.7 X 


1.4 


-63.4 


44.57 ± 0.03 


5.62 ± 


0.03 


4.05 ± 


0.07 


15.7 


3.6 


B 


D17.25-0.20b 




























D 


U18.15-0.28 


18.15 


-0.31 


2.2 


1.7 X 


0.7 


+77.9 


52.53 ± 0.04 


7.32 ± 


0.04 


4.75 ± 


0.06 


24.6 


7.6 


A 



Note. — Table[3]is published in its entirety in the electronic edition of the Astrophysical Journal Supplement Series. A portion is shown here for guidance 
regarding its form and content. 
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Table 4. Mean Properties of Nebulae with Associated ^^CO 





N 


\V\ Offset 


Tmb 


AV 


Size 


Tex 


Ar(i3CO) 






(kms-i) 


(K) 


(kms-l) 


(') 


(K) 


(xlO"* 


All 


253 


2.98 ± 2.41 


4.77 ± 2.32 


4.19 ± 1.42 


1.9 ± 1.3 


12.1 ± 4.8 


3.1 ± 2.6 


UC 


111 


3.03 ± 2.61 


5.21 ± 2.40 


4.43 ± 1.38 


1.7 ± 1.1 


12.2 ± 5.0 


3.5 ± 2.8 


Compact 


95 


3.01 ± 2.37 


4.96 ± 2.36 


4.23 ± 1.40 


2.2 ± 1.6 


13.3 ± 4.7 


3.3 ± 2.6 


Diffuse 


47 


2.81 ± 2.01 


3.32 ± 1.29 


3.56 ± 1.37 


1.9 ± 1.0 


9.3 ± 3.0 


1.6 ± 1.0 


A 


112 


2.65 ± 2.06 


5.63 ± 2.51 


4.52 ± 1.26 


1.7 ± 0.8 


12.8 ± 5.1 


4.0 ± 3.0 


B 


75 


3.00 ± 2.37 


4.62 ± 2.06 


4.05 ± 1.31 


2.1 ± 1.5 


12.1 ± 4.4 


2.8 ± 2.0 


C 


66 


3.53 ± 2.90 


3.46 ± 1.51 


3.80 ± 1.66 


2.2 ± 1.7 


11.0 ± 4.4 


1.9 ± 1.9 
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Fig. 1. — Longitude-LSR velocity diagram for Hll regions located inside the GRS survey 
zone. The nebulae are shown projected onto the Galactic plane. The symbols represent UC 
nebulae (small filled circles), compact nebulae (medium filled circles), and diffuse nebulae 
(large open circles). 
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Fig. 2. — The CO/Hll association procedure for the U43.24-0.05 Hll region (see text). 
Panel A: Using single velocity channels near the RRL velocity and nominal Hll region 
position we first locate the associated molecular emission. The black cross shows the Hll 
region position. The black rectangle defines the extent of the ^^CO emission that is deemed 
to be associated with the H ll region. Panel B: The average ^^CO spectrum of the GRS data 
cube voxels that lie within the panel A black rectange is shown. The solid vertical line fiags 
the H II region RRL velocity and the dashed lines show the velocity range of the associated 
emission. Panel C: The integrated intensity image made from the data cube using the line 
center velocity and FWHM line width found by the Panel B Gaussian fit is shown. This 
image is then used to define the final extent of ^^CO emission that is deemed to be associated 
with the Hll region. The black outhne shows this region which is defined by a threshold 
algorithm (see text). Panel D: The average ^^CO spectrum of the data cube voxels that 
lie within the panel C threshold defined region is shown. The vertical line flags the RRL 
velocity. The Gaussian fit to the emission line is used to derive the physical properties of 
the ^^CO gas associated with this Hll region. The spectrum from panel B is shown in gray 
for comparison. 
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Fig. 3. — Integrated intensity images (grey scale and black contours) for a sample of Hll 
regions with A confidence parameters. Nominal H ll region positions are marked with crosses 
and our fitted ellipses are also shown. Contours of 21cm VGPS continuum emission are shown 
in grey. Tickmarks on the VGPS contours point downhill, towards decreasing values. See §3 
for a detailed description of these images. 
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Fig. 4. — Integrated intensity images for a sample of H ll regions with B confidence param- 
eters. Other parameters are as in Figure O 
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Fig. 5. — Integrated intensity images for a sample of H ll regions with C confidence param- 
eters (top row), D (lower left) and E (lower right). Fitted ellipses are shown for C sources 
only. Other parameters are as in Figure [3l 
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Confidence Parameter 



Fig. 6. — The number distributions of the quahtative confidence parameter, CP. Sources 
with associated molecular gas are rated A, B or C in order of decreasing confidence in the 
association (see ^text). Confidence parameter D and E sources have little or no associated 
molecular emission, respectively. The top left panel is a stacked histogram where the top 
line shows the histogram for the entire sample of 301 Hll regions. The open, hatched and 
grey histograms show the contribution that UC, Compact, and Diffuse nebulae, respectively, 
make to the total in each bin. 
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Fig. 7. — Longitude-LSR velocity diagram for 253 Hll regions in our sample that show 
associated ^^CO emission. Symbols indicate the confidence parameter of each nebula: CP 
A (large filled circles), B (small filled circles), or C (small open circles). The lines show 
the locii of the LSR ter minal velocity ex pected from t wo different Galactic rotation curve 
models: Clemens Jl985h (dotted line) andlBrand] Jl986h (full line). 
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Tex (K) 



Fig. 8. — N e bular excitation temperature, Tex, derived from the ^^CO survey of 
Sanders et al. I (Il986l ) . Panels are the same as in Figure [61 
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Fig. 9. — Estimated H2 column density (top axis) and ^^CO column density (bottom axis) 
for the molecular clumps associated with our H ll regions. Panels are the same as in Figure 

la 
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Fig. 10. — The LSR velocity difference between the Hll region RRL and molecular gas ^^CO 
velocity. The distribution peaks near zero, as expected. Panels are the same as in Figure [6l 
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"•"mb (K) 

Fig. 11. — Main beam brightness temperature ^^CO line intensity from Gaussian fits. Panels 
are the same as in Figure El 
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Size (arcmin.) 



Fig. 13. — Angular size defined as the geometric mean diameter of tlie fitted ellipses. Panels 
are the same as in Figure El 
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Fig. 14. — ^'^CO line intensity plotted as a function of FWHM line width. Shown are 
all the CP A, B, or C sources: UC (filled triangles), Compact (filled circles), and Diffuse 
(open squares) nebulae. The lines dividing the figure into quadrants are from an analysis of 
molecular clumps in the GRS dataset (see §5.2p . 



